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Abstract 

Some time ago, Atiyah and Manton observed that computing the holonomy of 
Yang-Mills instantons yields good approximations to static Skyrmion solutions of the 
Skyrme model. This paper provides an extension and explanation of this result, by 
proving that instanton holonomies produce exact solutions of a BPS Skyrme model, 
in which the Skyrme field is coupled to a tower of vector mesons. Neglecting any 
(or indeed all) of the vector mesons breaks the scale invariance and removes the BPS 
property of the Skyrmions. However, it is shown that a truncation of the BPS Skyrme 
theory, in which only the first vector meson is included, already moves the Skyrme 
model significantly closer to the BPS system. A theory that is close to a BPS system 
is required to reproduce the experimental data on binding energies of nuclei. A zero- 
mode quantization of the Skyrmion is performed in the truncated BPS theory and the 
results are compared to the physical properties of the nucleon. The approach is an 
analogue in five-dimensional Minkowski spacetime of a recent holographic construction 
of a Skyrme model by Sakai and Sugimoto, based on a string theory derivation of a 
Yang-Mills-Chern-Simons theory in a curved five-dimensional spacetime. 



1 



1 Introduction 



The Skyrme model [1] is a nonlinear theory of pions with topological soliton solutions, 
called Skyrmions, that are identified as baryons. In this paper, three outstanding Skyrmion 
issues are drawn together and progress is made by introducing an extended Skyrme model, 
obtained from five-dimensional Yang-Mills theory. The three outstanding issues that the 
new extended theory addresses £1X6 clS follows. 

Firstly, the Skyrme model is not a BPS (Bogomolny-Prasad-Sommerfield) theory, in the 
sense that the soliton solutions do not attain the topological lower bound on the energy. 
In fact, the single Skyrmion exceeds the topological energy bound by 23% in the case of 
massless pions. This energy excess allows the possibility of a significant classical binding 
energy for higher charge Skyrmions, and indeed this is the case [2] : for example, the energy 
of the baryon number two Skyrmion exceeds the topological bound by 18%, which is already 
5% lower than the single Skyrmion. Such binding energies are much greater than those 
observed experimentally in nuclei, where binding energies are typically less than 1%. A 
BPS Skyrme model would therefore appear to be a better starting point for obtaining more 
realistic binding energies, since a small perturbation away from a BPS theory is likely to 
produce the required small binding energies. Motivated by this application a BPS Skyrme 
model is introduced, in which the usual Skyrme model is extended by the inclusion of an 
infinite tower of vector mesons. Furthermore, neglecting some of the vector mesons provides 
a natural way to perturb away from the BPS system. 

Secondly, Atiyah and Manton [3] have shown that computing the holonomy of Yang-Mills 
instantons yields good approximations to static Skyrmion solutions of the Skyrme model. It 
is not obvious why this approximation turns out to be so successful, or if instanton holonomies 
give exact solutions of some modified Skyrme model. This paper provides an explanation 
of the Atiyah-Manton procedure, by proving that the holonomy of Yang-Mills instantons 
yields exact solutions of the BPS Skyrme model. The BPS Skyrme model reverts to the 
usual Skyrme model if all the vector mesons are neglected, and this explains the accuracy 
of the Atiyah-Manton approximation. Furthermore, it is shown that including only the first 
vector meson already significantly improves the accuracy of the instanton approximation 
of the Skyrme field and allows an excellent approximation to the vector meson field to be 
extracted from the instanton. 

Thirdly, the usual Skyrme model includes pion degrees of freedom, but neglects all the 
other mesons. There is a long history of attempts to include other mesons, particularly the p 
meson [4, 5], but there are difficulties because of the large number of coupling constants that 
need to be determined: although some progress has been made using the ideas of hidden local 
symmetry and vector meson dominance [6, 7, 8]. In the BPS Skyrme model, all parameters 
are uniquely determined once the energy and length units are fixed. Truncating the BPS 
theory, for example by including only the first vector meson, allows the usual Skyrme model 
to be extended without the introduction of further unknown parameters. The Skyrmion is 
studied in this truncated BPS theory, including its zero-mode quantization, and the results 
are compared with the physical properties of the nucleon. 

The techniques used in this paper, to derive the BPS Skyrme model and its connection 
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to instanton holonomies, are inspired by the work of Sakai and Sugimoto [9]. Using a 
string theory construction and holographic methods, they were able to derive a Skyrme 
model coupled to an infinite tower of massive vector mesons from a Yang-Mills-Chern-Simons 
theory in a curved five- dimensional spacetime. The Skyrme field of their extended Skyrme 
model corresponds to the holonomy of the curved space instanton, though unfortunately 
this instanton solution has not yet been determined, even numerically. The Skyrmion in 
the truncated version of the extended Skyrme model, which includes only the pion and p 
meson degrees of freedom, has been investigated [10], though its quantization has not. A 
collective coordinate quantization of the instanton has been performed [11, 12], but only by 
approximating the true curved space Yang-Mills-Chern-Simons instanton by the flat space 
Yang-Mills instanton. 

In some respects, the work in the present paper may be regarded as a five- dimensional 
Minkowski spacetime analogue of the five- dimensional curved spacetime theory of Sakai and 
Sugimoto [9], with the advantage that the instanton, and various other ingredients, can 
be found explicitly. Of course, a disadvantage of this work is that there is no AdS/CFT 
correspondence to justify the approach, even though several results are qualitatively similar 
to those of the Sakai-Sugimoto theory, suggesting that there are some merits in considering 
this theory. 

2 Skyrmions and instantons 

In the Skyrme model [1] the pion degrees of freedom are encoded into an S77(2)-valued 
Skyrme field U. In the massless pion approximation, the static energy of the Skyrme model 
is 

£s ky = / ( - ^TriRiRi) - ^ Tr ([^ ^l 2 )) (2-1) 

where Ri = diU ' U~ x is the su(2) — valued current. In the above, e is the dimensionless 
Skyrme parameter and f w may be interpreted as the pion decay constant. Note that there 
are differing conventions, related by a factor of 2, for the pion decay constant. In this paper 
the convention is chosen so that the physical value is /„- = 92.6 MeV, which agrees with the 
convention in [9] and related papers. 

The parameters f n and e, whose values are to be fixed by comparison with experimental 
data, merely set the energy and length units and can be scaled away. Explicitly, if energy 
units of f n /2e and length units of l/ef w are used, then in dimensionless Skyrme units the 
energy becomes 

E Sky = J ( - ^Tr(i^) - ^Tr([i4 d 3 x. (2.2) 

This dimensionless form is used in the remainder of this section. 

The Skyrme field is required to tend to a constant element of SU(2) at spatial infin- 
ity (usually chosen to be the identity matrix) and this compactifies space to S* 3 . A given 
Skyrme field therefore has an associated integer topological charge B e Z = 7r 3 (577(2)) 
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given explicitly by 



1 




(2.3) 



B = 



24tt 2 



x. 



It is this topological charge that is to be identified with baryon number [13]. The Skyrmion 
of charge B is the field U that is the global minimum of the energy (2.2) for all fields in the 
given topological charge sector. 

The Faddeev-Bogomolny bound [14] states that 



and it is easy to prove that this bound cannot be attained for non-zero B. 

Recall that BPS solitons may be defined as solutions in which a topological energy bound 
is saturated and therefore, in this sense, Skyrmions are not BPS solitons. Skyrmion solutions 
can only be obtained numerically and, as mentioned in the previous section, the energy of 
the B = 1 Skyrmion is 127T 2 x 1.23 and the energy of the B = 2 Skyrmion is 247T 2 x 1.18. 
Numerical Skyrmion solutions have been obtained up to reasonably large baryon numbers 
[15] and reveal that the energies for larger values of B are significantly closer to the bound 
than for these low charge Skyrmions: for example, the B — 17 Skyrmion has an energy less 
than 17 x 127r 2 x 1.08. Computations based on periodic Skyrme fields [16] predict the limiting 
value E Sky /B — > 12-n 2 x 1.036, as B — > oo. The fact that the energy of the single Skyrmion 
is much further from the bound than for larger values of B implies binding energies that are 
much greater than those found experimentally for nuclei, which are typically less than 1%. 
The small binding energies of nuclei therefore motivate the search for a BPS Skyrme model, 
in which binding energies would vanish, allowing the possibility that a small perturbation of 
the BPS system might result in realistic nuclear binding energies. 

Although Skyrmions can only be obtained numerically, there are two analytic methods 
that produce Skyrme fields which are excellent approximations to the true Skyrmion solu- 
tions. One approach is the rational map approximation [17], but the method of interest in 
this paper is that of Atiyah and Manton [3] in which a Skyrme field is generated from the 
holonomy of a Yang-Mills instanton in R 4 . This is briefly reviewed below. 

Let Ai be the components of an SU(2) Yang-Mills instanton in M 4 , where uppercase latin 
indices run over all four space coordinates I = 1,2,3,4. The Skyrme field is defined to be 
the holonomy of this instanton computed along lines parallel to the X4-axis. Explicitly, 



where V denotes path ordering and x = (xi,x 2 ,x 3 ) are the Cartesian coordinates in the 
remaining M 3 C M 4 . As A 4 takes values in the Lie algebra su(2) its exponential is group- 
valued, so that C/(x) : R 3 >->■ SU(2), as required for a static Skyrme field. The ± factor in 
(2.5) is because the holonomy should really be defined on a closed loop on S* 4 and the sign 
may be required to account for the transition function that connects — oo to oo, corresponding 
to the same point on S* 4 . 



£ Sky > 12n 2 \B 



(2.4) 




(2.5) 
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As shown by Atiyah and Manton [3], the baryon number of this Skyrme field is equal to 
the instanton number of the gauge field, that is, B = N where 

N = J Tr(F u *F u )d 4 x, (2.6) 

and the dual field strength is defined by *Fjj = \eijklFkl- 

The Yang-Mills theory is conformally invariant and hence the instanton field includes 
an arbitrary scale. This construction does not provide an exact solution of the Skyrme 
model for any instanton, but for each N a suitable choice of instanton, including its scale, 
provides a remarkably good approximation to the static Skyrmion with baryon number N. 
The energy of the approximate Skyrme field is typically around a percent higher than that of 
the numerical Skyrmion and correctly reproduces the symmetry of the Skyrmion for a range 
of highly symmetric cases studied to date. For example, instantons have been constructed 
that correspond to the B = 1 spherically symmetric and B = 2 axially symmetric Skyrmions 
[3], tetrahedral and cubic Skyrmions with B = 3 and B = 4 [18], icosahedrally symmetric 
Skyrmions with B = 7 and B — 17 [19, 20], and the triply periodic Skyrme crystal [21]. 
There is therefore significant evidence to support the correspondence between Skyrmions and 
instanton holonomies, though a deeper understanding of this connection and its remarkable 
accuracy is desirable. 

Recently, the representation of a Skyrme field as an instanton holonomy has reappeared 
in the context of five- dimensional theories in compactified and/or curved spacetimes [22, 
23, 9, 24]. Although these approaches certainly have the flavour of the Atiyah-Manton 
construction, none of them actually involve exact self-dual Yang-Mills instantons in R 4 . 
However, as the Sakai-Sugimoto construction [9] is a main motivation for the present paper, 
and the techniques used here have some similarities to that work, it is perhaps useful to give 
a brief overview of the relevant aspects of this model. 

The Sakai-Sugimoto theory is based on a holographic approach to QCD in the limit of a 
large number of colours, using the AdS/CFT correspondence to map to a dual string theory 
consisting of probe D8-branes in a background of D4-branes. The action on the probe D8- 
branes leads to a Yang-Mills-Chern-Simons theory in a five- dimensional curved spacetime. 
The spacetime involves a warped product of (3+l)-dimensional Minkowski spacetime and 
an additional holographic direction. The static soliton in this theory is interpreted as the 
baryon, and has a fixed size determined by the ratio of the Chern-Simons coefficient to 
the curvature associated with the holographic direction. Unfortunately, the soliton in this 
theory has not been determined, even numerically. The expectation is that for a sufficiently 
small Chern-Simons coupling the soliton size will be small enough that the soliton can be 
approximated by a soliton of the flat space theory, which is simply a self-dual Yang-Mills 
instanton in K 4 , with a particular small scale. In fact, all work to date on this theory has 
used the flat space Yang-Mills instanton approximation, see for example [11, 12]. The dual 
theory in (3+l)-dimensional Minkowski spacetime is obtained by performing an expansion 
of the five-dimensional theory in terms of Kaluza-Klein modes of the holographic direction. 
These modes contain the Skyrme field plus an infinite tower of massive vector mesons and 
the associated theory is an extended Skyrme model, which reverts to the usual Skyrme 
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model if the massive vector mesons are ignored. In summary, the Sakai-Sugimoto theory 
provides a correspondence between Yang-Mills-Chern-Simons instantons on a curved four- 
manifold and an extended Skyrme model. One of the results of the present paper is to 
produce an analogous correspondence where the curved four-manifold is replaced by M 4 and 
the Chern-Simons term is absent. This allows a direct connection to be made to the Atiyah- 
Manton construction and also provides a natural extension of this method, together with an 
understanding of the remarkable accuracy of this approach. 

It should be noted that a realization of the Atiyah-Manton construction has been pro- 
posed [25] in which instantons appear as domain wall Skyrmions in a five- dimensional Yang- 
Mills-Higgs theory. This is certainly different from the approach discussed in the present 
paper and does not involve vector mesons. It also implies that terms involving higher deriva- 
tives than the Skyrme term should appear, which is not the case here, though it might be 
possible that some connection could be made between the two approaches by integrating out 
the vector mesons. 



As the details of the derivation for the full non-abelian gauge theory are quite cumbersome, 
it is useful to first consider a similar approach in a prototype abelian gauge theory, where 
the formulae are more manageable. 

Consider an abelian gauge theory in five-dimensional Minkowski spacetime with coordi- 
nates t, xi, where I = 1,2,3,4. For notational convenience define z = x<± and let lowercase 
latin indices run over the three remaining spatial coordinates % = 1,2,3. The real-valued 
gauge potential has components As most of this paper will be concerned with 

static fields with a t — 0, attention may be restricted to the static Yang-Mills energy 



where f LJ = <9/aj - <9ja/. 

As mentioned in the previous section, a crucial ingredient of the Sakai-Sugimoto con- 
struction [9] is the expansion of the gauge potential in terms of Kaluza-Klein modes in 
the holographic direction. In flat Euclidean space a replacement needs to be found for the 
curved space Kaluza-Klein modes. Simply taking the zero curvature limit is not suitable 
as the modes then degenerate to fourier modes, which are not appropriate on the infinite 
line. The required modes must form a complete orthonormal basis for square integrable 
functions on the real line with unit weight function (this is necessary to obtain canonical 
kinetic terms for the vector mesons). This problem is familiar to numerical analysts using 
spectral methods [26] and the recognized solution is provided by Hermite functions ip n {z), 
where n is a non-negative integer and 



3 An abelian prototype 




(3.1) 




(3.2) 
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In a gauge in which aj — > as \z\ — > oo, the components of the gauge potential can be 
expanded in terms of Hermite functions as 

oo oo 

a 2 (x, z) = J2 a"(x)^n ai(x, z) = E /^(x)^)- (3-3) 

ra=0 ra=0 

Consider a gauge transformation a/ >->■ 5/ = aj — djh, for which a z = 0. Clearly, this requires 
that d z h = a z , and hence h is given by 

J-oo n=Q \ J-oo J 

Hermite functions satisfy 



= \/f ^n-l(^) - y ^yV,+l(*), (3.5) 

where prime denotes differentiation with respect to z. This implies that 

f V2 P+ i(0^ = E^iK(4 (3.6) 

m=0 

f iMOdf = 7 2 ^ + W + E7 2 m p K + iW, (3.7) 

J-OO n 



m=0 



where 7^ and 7™ are non-zero constants. 

The additional function ip+{z) has been introduced and is defined by 

iM*) = "7^ f iMO # = J + \^{z/V2) (3.8) 

y27T4 J-00 z z 

with erf(z) the usual error function 

2 

0F 



erf(z) = — [ Z e-?dC (3.9) 

\/7T JO 



The normalization of ip+{z) has been chosen so that 00) = and vp + (oc) = 1. 

The gauge transformation (3.4) can now be written in terms of the basis functions 

ip + (z),i/j n (z) as 



fc(x, *) = «(X) + E ^( X ) ^n(^)- (3.10) 

n=0 

As V'n(oo) = and — 1? then -u(x) is identified as the holonomy 

/oo 
a,z(x, £) <i£. (3-11) 
-oo 
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In the new gauge, where a z = 0, then 

oo 

a t = a % - dih = -c^(x) ^j + {z) + £(/3?(x) - dih n (x)) ^ n (z). (3.12) 

n=0 

After defining the vector fields (x) = /3™(x) — <9j/i n (x) this becomes 

oo 

a % = -9 iU (x)^+(z) + X^<(x) Vn(^)- (3.13) 

n=0 

In this gauge the holonomy appears in the boundary condition Sj — >■ — diU as z — >■ oo. 
Using (3.13) and (3.5) the components of the field strength are 

oo 

j z . = _ d . u ^ +{z) + J2v^' n (z) (3.14) 



n=0 



and 

oo 

^^(9,<-^)^(z). (3.15) 



3 

n=0 



Using the orthonormality of the Hermite functions 

if m (z)ip n (z) dz = 5 mn (3.16) 
to perform the integration over z, the abelian Yang-Mills energy (3.1) becomes 



j — i 



E = J (^(^) 2 " ^i*" + Ej^W - ^<) 2 + ^n«) 2 " ^n<< +2 }) ^ 

(3-17) 

where the coefficients are m 2 n — n + \ and q n = ^ {n + l)(n + 2). 

The second term in (3.17) may seem a little strange, but it is simply the analogue in the 
prototype abelian theory of the familiar mixing between the Skyrme field and the lightest 
axial vector meson that arises in coupling the Skyrme model to vector mesons [5]. 

This approach has produced a correspondence between abelian Yang-Mills theory in R 4 
and a field theory in IR 3 containing an infinite tower of vector mesons, plus a scalar field 
related to the holonomy of the gauge potential. Note that m n are not the meson masses 
because q n ^ 0, hence the associated mass matrix is not diagonal. A truncated theory can 
be defined by including only the first Af vector mesons and setting i>™ = for all n > Af. 
The mass matrix is then diagonalized by an SO(Af) rotation of the remaining vector mesons 
and the meson masses determined from the eigenvalues of the Af x Af mass matrix. In the 
extreme case, Af — 1, where only the first vector meson remains, no rotation is required and 
the mass of this meson is obviously m = 1/ \/2. 
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Perhaps it is worth making a comparison between the above approach and the more 
common techniques of holographic QCD. In holographic QCD the curvature of the extra 
dimension induces a discrete spectrum and fields are then expanded in terms of the associated 
Kaluza-Klein modes. In the current situation the extra dimension is flat and therefore 
the spectrum is continuous. A discrete spectrum must be identified in order to mimic 
the holographic construction. The traditional approach to this problem in flat space is to 
compactify the extra dimension to produce a discrete spectrum. The continuous spectrum is 
then recovered in the limit of decompactification. However, for the application in the present 
paper it is not appropriate to compactify the extra dimension, because a modification of space 
means that the connection to the instanton in M 4 is then lost. Furthermore, the identification 
of the holonomy with the non-normalizable mode is no longer obvious in a compact extra 
dimension. The above Hermite truncation selects a discrete spectrum without the need to 
modify spacetime, and the continuous spectrum is recovered in the limit M — > oo. In this 
respect it is a crucial feature that the associated mass matrix is not diagonal. Increasing 
Af not only adds an additional mode but also shifts the masses of all the previous modes, 
allowing the continuous spectrum to reappear as Af — > oo. 

In principle, the eigenfunctions of any Sturm-Lioville operator on the line might be used 
as basis functions. However, a significant restriction is imposed by the requirement that the 
basis functions are orthonormal with respect to the unit weight function, which is needed 
to obtain the standard kinetic terms for the vector mesons. Hermite functions are a canon- 
ical choice and possess additional desirable features, such as the associated Sturm-Lioville 
operator involving only second and not first order derivatives. The precise combination of 
conditions that need to be imposed to uniquely arrive at the Hermite functions has not been 
investigated. 

In the next section a similar approach will be applied to SU (2) Yang-Mills theory, where 
it is shown that the S77(2)-valued scalar field that arises from the instanton holonomy is the 
Skyrme field, with energy function precisely that of the Skyrme model. 

4 The Skyrme model from Yang-Mills theory 

Consider SU(2) Yang-Mills gauge theory in R 4 with the sw(2)-valued gauge potential Aj, 
with 1 — 1,2, 3, 4, and again for convenience set z = x 4 . The Yang-Mills energy is given by 



where the factor of | is due to the normalization of the su(2) generators as — Tr(T a T&) = 25 a b- 
Starting with a gauge in which Aj — > as \z\ — > oo, the gauge A z = is obtained by 
applying the gauge transformation 




(4.1) 



Ai I— >■ gAig 1 -digg 



i 




with 




(4.3) 
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The holonomy is 

/oo 
Aj(x, z)dz, (4.4) 
-oo 

and in the new gauge the holonomy appears in the boundary condition for Ai, since now 
Ai — > —diU U' 1 as z — > oo. 

As in the abelian case (3.13), the gauge field can be expanded in terms of Hermite 
functions as 

oo 

JU = -diUU'^+iz) + W?( x ) M*)- ( 4 -5) 

n=0 

The emergence of the Skyrme model can be seen by first neglecting the vector fields W" '. 
With this truncation the components of the field strength are 

F zl = -diUU- 1 ^ = -Ri~r, Fij = [Ri, Rj]ijj + (ijj + - 1). (4.6) 

V27T4 



Substituting these expressions into the Yang-Mills energy (4.1), and performing the integra- 
tion over z, yields the energy of the Skyrme model 

E s = J ( - ^Tr(R t R t ) - ^Tr([R t , i^] 2 )) d 3 x, (4.7) 



where 

1 



ci = -r-f=, c 2 



/oo 
2ip 2 + {ip + — I) 2 dz — 0.198. (4.8) 
-oo 



This is the Skyrme model in dimensionless units, but it is not in Skyrme units because the 
constants c\ and c 2 are not unity. In these units the Faddeev-Bogomolny energy bound (2.4) 
becomes 

E s > 12tt 2 v ^c^ \B\ = 2.005 ir 2 \B\ . (4.9) 

This bound should be compared with the energy bound that derives from the full Yang-Mills 
theory, namely 

E>2tt 2 \B\, (4.10) 

which is attained by the instanton solutions. This shows that the two bounds are remarkably 
close, but that the Faddeev-Bogomolny bound is stricter by |%. Of course, the Faddeev- 
Bogomolny bound only applies to the usual Skyrme model, whereas the bound (4.10) is 
equally valid if some, or indeed all, of the vector mesons are included. 

The Skyrme model is not scale invariant, in contrast to the Yang-Mills theory, and hence 
the scale of the Skyrme model has emerged because of the truncation that ignores the vector 
mesons. Perhaps it is useful to think of this in terms of the theory flowing to a conformal 
theory as all the vector mesons are included. The issue of Skyrmion and instanton scales 
also appears in the Sakai-Sugimoto derivation [9] of the Skyrme model. In that case the 
Yang-Mills theory in curved space with a Chern-Simons term is not scale invariant and the 
instanton has a fixed small size, of the order of the string scale. The size of the Skyrmion 
in the usual Skyrme model is not related to the size of the instanton, but as more vector 
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mesons are included the size of the Skyrmion in the extended Skyrme model must tend to 
the small size of the instanton. 

The energy of a single Skyrmion exceeds the Faddeev-Bogomolny bound by about 23% 
and the energy of a Skyrme field generated from the holonomy of a single instanton exceeds 
the Faddeev-Bogomolny bound by about 24%, for an optimal choice of the instanton scale. 
When combined with the above bounds, this reveals that, for the optimal instanton scale, 
the result of neglecting all the vector mesons is to raise the energy by less than 25%. The 
vector meson terms in the expansion (4.5) have trivial topology and therefore the holonomy 
term captures all the topological features of the instanton (and hence the Skyrmion), but 
the above results demonstrate that it also captures most of the energetic properties too. 

Including the infinite tower of vector mesons produces a BPS Skyrme model, since the 
model is simply equivalent to Yang-Mills theory with one extra dimension. An infinite 
sequence of extended Skyrme models exist that interpolate between the usual Skyrme model 
and the BPS Skyrme model, as the number of included vector mesons ranges from zero to 
infinity. The remainder of this paper is devoted to a detailed analysis of the first member of 
this sequence that extends the usual Skyrme model. 



5 Including the first vector meson 

This section considers the extension of the Skyrme model obtained by including only the 
first vector meson, which physically corresponds to coupling the pion field to the p meson. 

The expansion (4.5) is not convenient once the vector mesons are included, because the 
fields do not have a definite parity. It is first necessary to perform an additional gauge 
transformation to obtain an expansion in terms of parity eigenstates. 

Given the holonomy U, define the SU(2)-vahied field S to be its positive square root, so 
that S 2 = U. After a gauge transformation by g — the expansion (4.5) takes the form 

1 1 oo 

A t = --P^ + -Q t + ]T V^ n (5.1) 

Z Z n=0 

where ip* = eri{z/y/2) and 

Pi = S- l diS + diSS-\ Q i = S- 1 d i S-d i SS- 1 . (5.2) 

The vector meson V™ is simply the previous vector meson Wp in the new gauge. 

Including only the first vector meson V t ° (and dropping the superscript for notational 
convenience) gives 



and 



F l3 = -[^^-(l-^ + ^-^o+^^o 2 

-([P^Vj] + [V h Pj])^ + ([Q i: Vj\ + [V t ,Q 3 ])^ . (5.4) 
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Substituting these expressions into the Yang-Mills energy and integrating over z produces 
the energy 

E = E s + Ey + E h (5.5) 
where E$ is the earlier Skyrme energy (4.7) and Ey is the vector meson energy 

Ey = J -Trj^V} - djVi) 2 + \m 2 V 2 + csidiVj - d^V,, Vj\ + c 4 [V h Vjf } d 3 x, (5.6) 

with mass m = and constants 

f°° 1 i f°° 1 1/1 

c 3=/ 7^0 dz = r , c 4 = / e*l>odz = -\—. (5.7) 



'-oo4 2^ 

In most phenomenological approaches to including the p meson, the energy Ey is taken 
to be that of a massive Yang-Mills field. However, Ey only has this form if 2C3/C4 is equal to 
unity, whereas 2c 2 / C4 = 2\/2/3 = 0.94, and hence there is a slight difference from a massive 
Yang-Mills theory. 

The interaction energy Ej is 

E l = l-Trj-cstP^K^V--^^ 

+^[Qi, ^](^ - d V t ) + c 3 [g,, ^][^, Vj\ + c 7 ([P 4 , ^] + [Vi, P 3 }f 
+^([Qi,V j \ + \y i ,Q s ]) 2 }cPx, (5.8) 



where the constants are 



/oo I 7j-'| /-oo ( 

^-(1-^)^0^=77^, c 6 = J_ oo -(l-^ 2 Mdz = 0.049 

r°° 1 11 
° 7 = 7-oc 32^ ^ = 32 " 2 ° 6 = °-°° 7 - (5 ' 9) 

The energy (5.5) corresponds to an extended Skyrme model in which both the pion and 
p meson degrees of freedom are included. Models similar to this have been investigated in 
the past [4, 5, 7, 8] but because of the difficulty in determining the many possible interaction 
coefficients some of these terms have not been included, and/or simple relations between 
the coefficients have been imposed, for example by generating the theory using principles 
of hidden local symmetry [6]. An advatange of the current derivation is that all interaction 
coefficients are uniquely fixed from the higher-dimensional Yang-Mills theory. 

A similar truncation of the Sakai-Sugimoto theory produces an energy of precisely the 
same form as that found here, but with different values for the coefficients, and the classical 
single Skyrmion has been studied numerically [10]. The results show that the classical energy 
of the Skyrmion in that case is reduced by about 10% in comparison to the usual Skyrmion, 
but no quantization of the Skyrmion has been performed. In a later section it will be shown 
that the energy of a single Skyrmion in the theory (5.5) is even lower than this, and moreover 
a quantization of the Skyrmion reveals that it is vital to take into account the quantum spin 
energy of the Skyrmion when attempting to match to the physical properties of the nucleon. 
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6 The Skyrmion and its instanton approximation 



As the extended model (5.5) is obtained directly from Yang-Mills theory then good ap- 
proximations to the energy minimizing fields of this model will be obtained by extracting the 
appropriate components of the instanton, in terms of the basis expansion used for the trun- 
cation. In this section this extraction is performed for the case of a single Skyrmion through 
a detailed calculation of the decomposition of the instanton in terms of the expansion (5.1). 

The single Skyrmion and the single instanton both have 5*0(3) symmetry, associated 
with spherical symmetry in IR 3 . It is therefore useful to introduce the following symmetric 
tensors 

where Xi = Xj/|x|. 

The N — 1 instanton located at the origin in M 4 is given by Ai = iA Ia r a where r a are 
the Pauli matrices and 

A ia = r](X ia + Y ia ) + (Z ia , A za = (Y ia , (6.2) 

with the functions rj and ( defined to be 

Z T 

X 2 + r 2 + z 2 X 2 + r 2 + z 2 

Here A is the arbitrary scale of the instanton. 

The gauge A z = is obtained after the gauge transformation 

g{x, z) = exp(iF(r, z)x a r a ), (6.4) 

where 

F(r, z) = £ dr, () « = -^{i + ^tan- (-^_,) }. (6.5) 
The associated holonomy has the standard hedgehog form 

U = exp(i/(r)x a r a ), (6.6) 

with profile function 

/(r) = F(r,oo) = - 7 ^-,. (6.7) 

The Skyrme field must tend to a constant element of SU(2) as r — > oo and this is usually 
taken to be the identity matrix, corresponding to the boundary condition /(oo) = 0. In 
this section it is slightly more convenient to take the non-standard choice that U — > — 1 as 
r — > oo, with the associated profile function boundary conditions /(0) = and /(oo) = — ir, 
as satisfied by the profile function in (6.7). This is not an important change but does mean 
that some additional factors of n do not need to be introduced and carried throughout the 
following calculation. 
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Figure 1: The Skyrme field profile function minimizing the energy Eg. The solid curve is the 
numerical solution and the dashed curve is the instanton approximation. 



First, consider restricting to the usual Skyrme energy E$, which will reproduce the results 
of Atiyah and Manton [3]. 

For a Skyrme field of the hedgehog form (6.6) the usual Skyrme energy (4.7) reduces to 
the expression 

A numerical minimization of E s yields an energy of E s = 1.236 x 2tt 2 , with the associated 
profile function f(r) displayed as the solid curve in Figure 1. 

Restricting to the instanton approximation (6.7), the usual Skyrme energy Es/2ir 2 is 
plotted as a function of the instanton scale A as the dashed curve in Figure 2. This com- 
putation reveals that within the instanton approximation E$ is minimized by an instanton 
with scale A = 1.72, which has an energy E s = 1.246 x 2ir 2 . For comparison, the instanton 
generated profile function with minimizing scale is displayed as the dashed curve in Figure 1. 
The instanton is found to generate a good approximation to the Skyrme field, with energy 
only 1% above that of the numerical solution. 

The next step is to extend these calculations to determine the vector meson fields hidden 
inside the instanton. The energy minimizing Skyrme field and vector meson field are then 
computed in the extended Skyrme model (5.5) and compared with those from the instanton 
approximation. 

Performing the gauge transformation (6.4) yields the gauge A z = 0, and in this gauge 
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Figure 2: The usual Skyrme energy E^j2n 2 (dashed curve) and the energy in the extended 
model E/27i 2 (solid curve) as a function of the instanton scale A. 



the remaining components are 

A ia = X ia cos 2F - (( + ±- ) sin 2f) + Y ia (77 -d r F) + Z ia (77 sin 2F + (( + ^-) cos 2F - . 

(6.9) 

Now perform the additional gauge transformation with 



g = S - l =exp(- t -fx a T a ), (6.10) 



to obtain 



A ia = X ia (r) cos H - (C + ^) sin if) + F ia (77 - ±d r H) + Z ia ^ sin H+(( + ^) cos if - , 

(6.11) 

where 

tf = ; 2r tan' 1 !" - = V (6.12) 

VA 2 + r 2 Vv/A 2 + r 2 / V ; 

The gauge potential (6.11) now has the correct parity properties to be compared with the 
expansion (5.1) in terms of parity eigenstates. An immediate comparison yields that 

1 1 

P ia = X ia -sinf + Yiaf , and Q ia = Z ia -(cos f - 1), (6.13) 

r r 

where Pi = iPi a T a etc. Then, by definition of the terms in the expansion (5.1), this gives 
that 

V£ = KX ia + k n 2 Y ia + k%Z ia , (6.14) 
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where the profile functions fc"(r) are given by the integrals 

fc"(r) = / [77 cos if - (C+ — ) sin if + — sin f i/jA ijj n dz (6.15) 

fcj(r) = \p-\d r H+^f^ n dz 

fc£(r) = ^j^sinii+^ + ^cosii-^cos/l^^. 

For the even vector mesons it is easy to show that k 2n (r) = /c| n (r) = due to the sym- 
metry ip2n(—z) = 4>2n(z), and for the odd vector mesons k 2n+1 (r) = due to the symmetry 
■02n+i( — z) = — ip2n+i(z). This is the correct parity associated with the fact that V™ is a 
vector meson for even n and an axial vector meson for odd n. 

For the first vector meson n — 0, and again dropping the superscript on V®, the above 
results reduce to 

V ia = p(r)Z ia , (6.16) 

where the profile function p(r) is 

p(r) — I irj sin H + (( + — ) cos H — --- cos fX ip dz. (6-17) 

The profile function p(r) satisfies the boundary conditions p(0) = p(oo) = 0. 

In terms of arbitrary profile functions f(r) and p(r) appearing in the spherical ansatz 
(6.6) and (6.16), the additional terms in the energy E = E s + Ey + Ej become 

Ey = 4tt r \p' 2 + % + ^ + m 2 p 2 + c 3 16^ + c 4 16p 4 ) r 2 dr, (6.18) 
Jo i. r r r J 

and 

El = 4,f { - 16c 5 ^(/ V + £) + £-5^) + ^(cos/ - 1) (6,9, 

-16c 6 ^-^ + 16c 3 ^(cos/ - 1) + 32c 7 /V + 2^(cos/ - l) 2 } r 2 

A numerical minimization of the extended energy E yields the value E = 1.060 x 2ir 2 . This 
shows that the energy of a Skyrmion in this theory is significantly closer to the topological 
lower energy bound (4.10) than in the usual Skyrme model. This result reveals that the 
truncation of the BPS Skyrme theory, in which only the first vector meson is included, already 
moves the usual Skyrme model significantly closer to the BPS theory. The numerically 
determined Skyrme profile function /(r) is displayed as the solid curve in Figure 3 and the 
vector meson profile function p(r) as the solid curve in Figure 4. 

Applying the instanton approximation, with profile functions (6.7) and (6.17), the energy 
E/2tt 2 is plotted, as a function of the instanton scale A, as the solid curve in Figure 2. The 
mimimizing instanton scale is A = 1.20 at which the energy is E = 1.071 x 2n 2 . This 
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Figure 3: The Skyrme field profile function minimizing the energy E. The solid curve is the 
numerical solution and the dashed curve is the instanton approximation. 



demonstrates that the instanton generated fields provide an excellent approximation to both 
the Skyrme field and the vector meson field, with an energy only 1% above that of the 
numerical solution. For comparison to the numerical solutions, the minimizing instanton 
profiles are displayed as the dashed curves in Figure 3 and Figure 4. Comparing Figure 1 
with Figure 3 confirms that, as expected, the instanton approximation to the Skyrme field 
profile function is much closer to the numerical solution in the extended Skyrme theory than 
in the usual Skyrme model: in fact, it is difficult to distinguish the two curves in Figure 3. 

7 Quantization of the Skyrmion 

It is well-known that fixing the parameters of the Skyrme model in the meson sector does 
not produce good results in the baryon sector. It is common practice to treat (at least some 
of) the meson constants as free parameters which are then fixed by comparison to selected 
baryon properties [27, 28, 4, 29, 30, 31]. It has been suggested [32] that the meson parameters 
could be interpreted as renormalized constants in the baryon sector, that result from known 
quantum effects not addressed within a simple zero-mode quantization. As demonstrated 
below, the same situation persists in the extended Skyrme model. 

Given the dimensionless formulation of the theory (5.5), the energy unit e and the length 
unit I are related to the pion decay constant f n and p meson mass m p by 

e = tfV2^/m p , I = l/(v / 2m p ). (7.20) 

This can be seen by performing the scaling E i->- eE, >->■ Xi/l, Vi >->■ ViJl/e, after which 
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Figure 4: The vector meson profile function p minimizing the energy E. The solid curve is 
the numerical solution and the dashed curve is the instanton approximation. 



the dimensionless Skyrme energy (4.7) takes the standard form (2.1) with Skyrme parameter 



m. 



e = (7.21) 



The vector meson energy (5.6) becomes 
Ey = J -Treaty - d.Vi) 2 + l -m 2 p V? + c 3 (<9^ - d 3 Vi)[V t , V 3 ] + c 4 [V t , V$] 2 } d 3 x, (7.22) 



where 



m~ 



c 3 = JL , c 4 = (7.23) 

For completeness, the additional coefficients in the interaction energy E\ remain unchanged 
except for the replacement 

The classical Skyrmion energy is E — eM, where M = 1.06 x 2-k 2 is the dimensionless 
static energy. Taking the physical values f n = 92.6 MeV and m p = 776 MeV, to set the units 
(7.20), gives a classical Skyrmion energy E = 580 MeV, which is far too low in comparison to 
the nucleon mass of 939 MeV. There is also a quantum contribution to the energy associated 
with the spin of the nucleon, but even if this contribution raised the total energy to that 
of the nucleon mass this would still not be an acceptable result, since physically the spin 
contribution needs to provide only a small contribution to the total energy. The quantum 
spin energy is calculated later in this section and reveals that taking the physical values for 
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the meson parameters f n and m p yields a quantum spin contribution of 1039 MeV, which is 
almost twice that of the classical energy and gives a total energy which is far too large. This 
shows that taking the physical values for the meson parameters does not produce acceptable 
results for the baryon. Confirmation of this is provided by calculating the size of the baryon, 
as follows. 

The physical value of the nucleon isoscalar root mean square radius is V < r 2 > = 0.72 fm. 
For the Skyrmion its dimensionless form R is calculated from the radial baryon density B as 

/•OO PCX) 9 

R 2 = r 2 Bdr = - r 2 -f sin / dr = 0.82. (7.25) 

JO JO 7T 

Inserting the length unit gives 

R 

y/ <r 2 > = — — . (7.26) 
V2m p 

Taking the physical value for m p , and using the fact that in natural units MeV" 1 = 197 fm, 
yields V < r 2 > = 0.16 fm, which is far too small. This is the origin of the excessive quantum 
spin energy mentioned above, when physical values are taken for the meson parameters. The 
baryon is far too small and hence so is its moment of inertia, which occurs in the denominator 
of the quantum spin energy. 

From now on the common practice is adopted of treating the meson parameters of the 
theory (in this case /„- and m p ) as free parameters that are to be fixed by comparison to 
physical properties in the baryon sector. 

The choice made in this paper is to fix the energy and length units by matching to the 
physical values of the nucleon mass and the isoscalar root mean square radius. Matching 
the latter, using equation (7.26), yields m p = 176 MeV, which is therefore only around a 
quarter of its physical value. To determine it is first necessary to calculate the quantum 
spin contribution to the nucleon mass, which is presented below. 

So far in this paper the discussion has been restricted to static fields. It is a simple matter 
to obtain the relevant Lagrangians from the static energies presented earlier by applying the 
obvious relativistic generalization. From the dimensionless form (5.5) of the static energy of 
the extended Skyrme model, the associated dimensionless kinetic energy is T — T s + T v + Tj 
where 

T s = /-Trjl^ + I^o,^] 2 }^, (7.27) 

T v = J -Tr|i(9 ^ - d.Vo) 2 + \m 2 V 2 + 2c 3 («9 ^ - d t V )[V , V,} + 2c 4 [V , d 3 x, 

(7.28) 

Ti = J -Trj - 2c 5 [P , PiMdoV - dM>) - 2c 6 [P , Pi\[V , V % ] - c 5 [P , PJ([Q , K] + Wo, Qi\) 
+\([Qo,Vi\ + Wo,Qi])(doVi - diV ) + c 3 ([Qo,^] + Wo,Qi])Wo,Vi\ 
+2c 7 ([P , Vi] + Wo, P^]) 2 + ^(lQo, V l ] + Wo, Qi\) 2 } d 3 x. (7.29) 
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The zero-mode quantization involves the rigid rotor ansatz 

S = e* SH Se-* tH , V^e^Vie'^, V = e^%e~^\ (7.30) 

where Q = iVL a r a is a constant element of su(2) determining the rotation frequency and 
axis, and S, Vi are the earlier static fields that minimize the classical static energy E. For 
time-dependent fields Vq can no longer be set to zero, as this is not consistent with the Gauss 
law for this system, that is, the field equation for Vq. The Gauss law requires that Vo has the 
same tensorial structure as [V, [V, Q]]. This determines the form of V = iV 0a T a to be [4] 

V 0a = xA + X2X a X b ^b, (7.31) 

where Xi( r ) i X2(j') are two additional radial profile functions with boundary conditions 
x'i(0) = Xi(°°) = X2(0) = ^2(00) = 0. The kinetic energy takes the form T = i/|f2| 2 , 
where \Q\ 2 = — |Tr(f2 2 ) and / = I s + J v + h is the moment of inertia, which after a tedious 
calculation is found to be 

/s = f /;{ Cl + C2 (r + ^)}sinV. 2 ^ (7.32) 



h = y j ( 00 {4p 2 + 4^ + 3x , 1 2 + 2 X / lX / 2 + X / 2 2 + m 2 (3x? + 2 XlX2 + X 2 ) 

+32c 3 p(2p X i + PX2 + + 64c 4 p 2 (2 X 2 + 2 X i X 2 + xl))r 2 dr, (7.33) 



T Jo I 32c s sm f[-^— (P ~ —) + fXiJ - 64cep— — xi 

o • 2/ ( q ^m 2 f , Xi\ . 2 , o Xi . x\ , Q /n Xi . n XiX2 ,x\, n\ 

-8 sin - - 8c 5 (p + — ) + p + 2p— + + 8c 3 p(2— + 2 + — + Xip) 

2 \ r r r r z rrr j 

1 RA ( f*2 2 , ■ 2f(( Xl\2 , (X1 + X2) 2 ' 

+64c 7 (/ Xi + sm fUp-—) + — 



sin 4 ^ Up+-) 2 + KX1 V ^ ^ ( 7 - 34 ) 



4/// , Xh2 , (X1 + X 2) 2 

2T r j r 2 

The functions X11X2 are determined by the V field equation and this is equivalent to 
the minimization of iy + h, given the profile functions / and p. The minimizing profile 
functions are presented in Figure 5 and the associated moment of inertia is computed to be 
I = l s + (/ v + ij) = 13.73 + 1.96 = 15.69 

In terms of the spin J = I\Q\ the dimensionless quantum spin energy is 

E Q = (7.35) 
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Figure 5: The profile functions Xi (solid curve) and X2 (dashed curve). 



where J 2 = j(j + 1) and the quantum spin number j — \ for the nucleon. The moment of 
inertia / has units of el 2 and the quantum spin energy has units which are the reciprocal of 
this. The mass of the nucleon is therefore 

1 f 2 m 3 l~2 3 

Taking the previously determined value m p = 176 MeV and requiring that (7.36) reproduces 
the nucleon mass M N = 939 MeV yields f n = 55.1 MeV, which is around 60% of its physical 
value. With these parameter values the classical and quantum spin contributions to the 
nucleon mass are 905 MeV and 34 MeV respectively, which is an acceptable split. 

By equation (7.21), this set of parameter values gives a Skyrme parameter e = 3.81. It is 
interesting to note that these parameters are reasonably close to those suggested in the usual 
Skyrme model by fitting to the properties of the lithium-6 nucleus, which give f n = 37.6 MeV 
and e = 3.26 [30]. 

8 Conclusion 

Inspired by methods of holographic QCD, a sequence of extended Skyrme models has been 
introduced that interpolate between the usual Skyrme model and a BPS Skyrme model. This 
provides an explanation and extension of the Atiyah-Manton construction of Skyrme fields 
from instanton holonomies, as this construction produces exact solutions of the BPS Skyrme 
model. 

The first extended Skyrme model is a nonlinear theory of pions coupled to the p meson and 
this has been investigated in some detail. The results reveal that this model is significantly 
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closer to a BPS theory than the usual Skyrme model, and it has been demonstrated that an 
extension of the Atiyah-Manton construction provides an excellent approximation to both 
the pion and p meson fields. This is encouraging as the experimental data on nuclear binding 
energies reveals that they are typically less than 1%, suggesting that a model close to a BPS 
theory is required. 

There are several avenues for future research that follow from the work presented in this 
paper. An obvious next step is to investigate multi-Skyrmions in the extended theory, to 
confirm and evaluate the reduced binding energies. This could be performed either using 
numerical full field simulations, similar to those applied in the usual Skyrme model [15], or 
using the instanton approximation. A study of the extended Skyrme models that include 
additional vector mesons would also be of interest: in particular it would be useful to compute 
the change in binding energies as more vector mesons are introduced. 

Many aspects of the Sakai-Sugimoto theory have been investigated since its introduction 
and it would be interesting to attempt similar studies for the flat space analogue introduced 
here. Examples of aspects to study include finite baryon density and temperature [33, 
34, 35, 36, 37, 38], the determination of interaction coupling constants and form factors 
[39, 40, 41, 42], together with an analysis of the nuclear force [43]. 

The Skyrme models considered in this paper are all applicable to massless pions, but 
it has been shown that in the usual Skyrme model there are significant differences in the 
massive pion theory, and the differences are encouraging in respect to comparisons with 
the properties of nuclei [44, 45, 31]. A pion mass term could simply be included in the 
extended Skyrme models, although the connection to a BPS Skyrme theory would then be 
lost. However, it is still to be expected that binding energies would be reduced in comparison 
to the usual Skyrme model, as the resulting increase in energy applies to Skyrmions of all 
baryon numbers. The Atiyah-Manton construction does not produce Skyrme fields with 
asymptotic fields appropriate to massive pions, but a modification of this construction has 
been introduced for massive pions [46], based on a connection to hyperbolic Skyrmions. It 
would be interesting to see if this modified construction can be understood in terms of the 
techniques introduced in the present paper, and perhaps this might lead to an extended BPS 
Skyrme theory for massive pions. 

Finally, during the preparation of this manuscript a preprint appeared [47] introducing 
novel Skyrmions in a different BPS Skyrme model. The BPS Skyrme model in question is 
of the type introduced some time ago [48], involving only a pion mass term and a term of 
sixth-order in the derivatives of the Skyrme field, and is therefore quite different from the 
one considered in the present paper. The novel Skyrmions presented in [47] are examples of 
compactons, that is, they have compact support. Such Skyrmions are trivially BPS solitons, 
since compactons placed far enough apart do not interact at all. However, there is an 
interesting mathematical structure, associated with an infinite dimensional symmetry, that 
for example allows the single Skyrmion to be obtained as a solution of a first order equation. 
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